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THE ROLE OF BAG SURFACE TENSION IN COLOR CONFINEMENT 

K. A. Bugaev^, A. I. Ivanytskyi^ 

Bogolyubov Institute for Theoretical Physics, Kiev, Ukraine 

We discuss here the novel view at the color confinement which, on the one hand, allows us to find out the 

surface tension coefficient of quark gluon bags and, under a plausible assumption, to determine the endpoint 

temperature of the QCD phase diagram, on the other hand. The present model considers the confining color 

^^ tube as the cylindrical quark gluon bag with non-zero surface tension. A close inspection of the free energies of 

^_^ elongated cylindrical bag and the confining color tube that connects the static quark-antiquark pair allows us 

^^ to find out the string tension in terms of the surface tension, thermal pressure and the bag radius. Using the 

C^ derived relation it is possible to estimate the bag surface tension at zero temperature directly from the lattice 

QCD data and to estimate the (tri) critical endpoint temperature. In the present analysis the topological free 

energy of the cylindrical bag is accounted for the first time. The requirement of positive entropy density of 

such bags leads to negative values of the surface tension coefficient of quark gluon bags at the cross-over region, 

i.e. at the continuous transition to deconfined quarks and gluons. We argue that the cross-over existence at 

I— I supercritical temperatures in ordinary liquids is also provided by the negative surface tension coefficient values. 

1-M It is shown that the confining tube model naturally accounts for an existence of a very pronounced surprising 

^~^ maximum of the tube entropy observed in the lattice QCD simulations, which, as we argue, signals about the 

Q^ fractional surface formation of the confining tube. In addition, using the developed formalism we suggest the 

1^ gas of free tubes model and demonstrate that it contains two phases. 



I— » 



1 Introduction 



^ A new paradigm of heavy ion phenomenology that the quark gluon plasma (QGP) is a strongly interacting 

^^ liquid [1 proved to be very successful not only in describing some of its properties measured by lattice quantum 

_!. chromodynamics (QCD), but also in explaining some experimental observables that cannot be reproduced 

^-H otherwise. Probably, the two most striking conclusions obtained within the new paradigm are as follows: first, 

1^ at the cross-over temperature, where the string tension of color tube is almost vanishing, the potential energy 

^^ of color charge is of the order of a few GeV [2 , i.e. it is 10 times larger than its kinetic energy, and, second, 

^-H the QGP, so far, is the most perfect fluid since its shear viscosity in units of the entropy density is found to 

^"H be the smallest one [H [4] . The first of these conclusions tells us that at the cross-over region there is no color 

^ charge separation [5] , whereas the second one naturally explains the great success of ideal hydrodynamics when 

'kJ( applied to relativist ic heavy ion collisions. 

rS Here we would like to discuss the recent progress achieved in our understanding of both the confinement 

^ phenomenon [6l [7] and the physical origin of the cross-over [HI |9l [TOl [TTJ [12] . As we demonstrate below such 

a progress was made possible after realizing a principal role played by negative values of the surface tension 
coefficient of large QGP bags pi [8]. Also here we argue that the negative values of the surface tension, that are 
responsible for an existence of the cross-over transition to QGP at low baryonic densities, play the same role in 
ordinary liquids. Moreover, in this work we would like to draw an attention to the problem of the temperature 
dependence of surface tension coefficient in liquids by clearly showing that for many liquids the well known 
Guggenheim relation (see Eq. ([To])) is not so well established experimentally as it is usually believed. 

The work is organized as follows. Section 2 is devoted to the confining tube model, in which the Fisher 
topological term of the QGP bag free energy is accounted for. In section 3 we show that at the cross-over region 
the surface tension of QGP bags is necessarily negative and argue that this is the case for ordinary liquids as 
well. The maximum of the tube entropy observed in lattice QCD is explained in section 4, where the model of 
gas of free tubes is also developed. The conclusions are given in the last section. 

2 Color confining tube and sQGP 

A color confinement, i.e. an absence of free color charges, is usually described by the free energy of heavy 
(static) quark-antiquark pair Fqq(T, L) = cTstr - L. In the lattice QCD the functional dependence of Fqq{T, L) 
on the temperature T and the separation distance L can be extracted from the Polyakov line correlation in a 
color singlet channel. Then it is customary to define 
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Figure 1. The internal energy Uoo (left) and entropy ^oo (right) of confining tube connecting two static color 
charges found by lattice QCD simulations for infinite separation distance between charges [2 . The internal 
energy is shown for 2 quark fiavors. 



• confinement: the case of non-zero string tension, i.e. (Tstr > 0; 

• deconfinement: the case of vanishing string tension astr -^ at T ^ Tco-, but one should remember 
that there is no color charge separation up to T > 1.3 Tco values of the cross-over temperature Tco [BIS]- 

The explanation of the latter is as follows: although at large distances L ^ oo the potential energy of static 



qq pair is finite Uqq{T^L) = Fqq — T-^jr- = Fqq -\- TSqq near Tco-, the values of /7(T, oo) are very large (see 
Fig. fl]). From Fig. IT] one can conclude that near Tco region QGP is a strongly interacting plasma (sQGP) which 
is similar to a liquid, since the ratio of the quark potential energy to its kinetic energy, the so called plasma 
parameter, ^^^ G 1 — 10 has the range of values that is typical for ordinary liquids [1 . 

The second striking feature of the confining tube can be seen in the right panel of Fig. [l] which clearly 
demonstrates that at T = Tco the entropy of static qq pair is very large SqqiTco^ oo) '^ 20. Such a value signals 
that really a huge number of degrees of freedom ~ exp(20) is involved, but the origin of large energy 
UqqiT^ oo) and entropy SqqiT^ oo) values near Tco foi" awhile remained mysterious [1 despite many attempts 
to explain it. 

Another problem of principal importance for phenomenological models of deconfinement phase transition 
[HI [9l [TOl [m [121 [El m [IS] is the value of the surface tension coefficient (J surf of QGP bags. There are several 
estimates for the surface tension coefficient (J surf of QGP bags [16 , but the question is whether can we determine 
cFsurf from lattice QCD? Therefore, in this section we consider an approach that allows us to determine the 
surface tension coefficient of QGP bags directly from the lattice QCD. As it will be shown in the section 4 such 
an approach naturally explains an existence of the 'mysterious maximum' [1 of the confining tube entropy. 

In order to estimate the surface tension of QGP bags let us consider the static quark-antiquark pair connected 
by the unbreakable color tube of length L and radius R <^ L. In the limit of large L the free energy of the color 
tube is Fstr -^ cFstrL. Now we consider the same tube as an elongated cylinder of the same radius and length 
[6]. In this case we neglect the free energy of the regions around the color charges, but for our treatment of large 
separation distances L ^ R this is sufficient. For the cylinder free energy we use the standard parameterization 
l8l[a[lQl[IIl[l2] 



Fcyi{T,L,R) = -p^{T)7rR^L ^ cTsurf{T)27rRL ^Tr\n 



ttR^L 



Vo 



(1) 



where Pv{T) is the bulk pressure inside a bag, (Jsurf{T) is the temperature dependent surface tension coefficient, 
while the last term on the right hand side above is the Fisher topological term [17 which is proportional to 
the Fisher exponent r = const > 1 [3 O [TOl [Til [121 [T5l and Vb ~ 1 fm^ is a normalization constant. Since we 
consider the same object then its free energies calculated as the color tube and as the cylindrical bag should be 
equal to each other. Then for large separating distances L ^ R one finds the following relation 



CFstr{T) = CFsurf iT) 271 R - Py{T)7TR^ 



Tr 



In 



ttR'^L^ 



(2) 



In doing so, in fact, we match an ensemble of all string shapes of fixed L to a mean elongated cylinder, which 
according to the original Fisher idea [17 and the results of the Hills and Dales Model (HDM) [18, 19 represents 
a sum of all surface deformations of such a bag. The last equation allows one to determine the T-dependence 
of bag surface tension as 



rf{T) 



CTstrjT) 1 

27ri? 2 



Pv{T)R 



Tt 
2ttRL 



In 



'ttR'^L 



(3) 
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if R{T)^ (Jstr{T) and Pv{T) are known. This relation opens a principal possibility to determine the bag surface 
tension directly from the lattice QCD simulations for any T. Also it allows us to estimate the surface tension 
at T = 0. Thus, taking the typical value of the bag model pressure which is used in hadronic spectroscopy 
Py{T = 0) = -(0.25)'^ GeV^ and inserting into ^ the lattice QCD values R = 0.5 fm and astr{T = 0) = (0.42)^ 
GeV^ [20 , one finds asurf{T = 0) = (0.2229 GeV)^ + O.bpyR ^ (0.183 GeV)^ ^ 157.4 MeV fm-^ [6 . The 
last term in Q does not modify our above estimate at T = 0, but, in contrast to [H [7], we keep it in order to 
demonstrate its importance for the confining tube with free color charges. 

The found value of the bag surface tension at zero temperature is very important for the phenomenological 
equations of state of strongly interacting matter in two respects. Firstly, according to HDM the obtained value 
defines the temperature at which the bag surface tension coefficient changes the sign [181 [HI El 

Ta = (Jsurf(T = 0) V} ■ X-^ e [148.4; 157.4] MeV , (4) 

where the constant A = 1 for the Fisher parameterization of the T-dependent surface tension coefficient ^7\ 
or A ^ 1.06009, if we use the parameterization derived within the HDM for surface deformations [18l [191 [Zj- 
Secondly, according to one of the most successful models of liquid-gas phase transition, i.e. the Fisher droplet 
model (FDM) [17 the surface tension coefficient linearly depends on temperature. This conclusion is well 
supported by HDM and by microscopic models of vapor-liquid interfaces [21 . Therefore, the temperature T^r 
in Q, at which the surface tension coefficient vanishes, is also the temperature of the (tri)critical endpoint Tcep 
of the liquid-gas phase diagram. On the basis of these arguments in Ref. [7 we concluded that the value of 
QCD critical endpoint temperature is Tcep = ^a = 152.9 ±4.5 MeV. Hopefully, the latter can be verified by the 
lattice QCD simulations using Eq. ([3|. 

Now the question is what is the surface tension coefficient above Tcep^ i.e. at supercritical temperatures. 
There are no experimental data on usual liquids in this region. In FDM and in the other well known model 
of liquid-gas phase transition, the statistical multifragmentation model (SMM) [22][23l[2T, the surface tension 
at supercritical temperatures is assumed to be zero, while in other models such a question is usually not 
discussed. The only exceptions known to us are the exactly solvable statistical models of quark gluon bags with 
surface tension [8^ 9 , their extension which includes the finite widths of large/heavy QGP bags [10, llf il2] 
and recently formulated generalization of the SMM [25 . For all these models it was demonstrated that the 
negative surface tension is the only physical reason of degeneration of the 1-st order phase transition into cross- 
over at supercritical temperatures. The question is whether the above suggested formalism can support such a 
conclusion. 

3 Surface tension coefficient at the cross-over temperature 

The above results, indeed, allow us to tune the interrelation with the color tube model and to study the bag 
surface tension near the cross-over to QGP. Consider first the vanishing baryonic densities. The lattice QCD 
data indicate that at large R the string tension behaves as 

CTstr = -^ ^0, (5) 

where Lq > and ^o > are some positive constants. Assuming the validity of Eq. ^ in the infinite available 
volume, one finds that for cFstriT) -^ +0 the string radius diverges, i.e. R ^ oo. 

Using Eqs. dl]) and O) we can write the total pressure ptot of the cylinder as follows 



pUL,R,T) = MT)-"-^^^^ ^^ --surfiT) a.,r , Tr 



R ttR^L R ttR^ ttR^L 

f{T) goln{L/Lo) , Tt 



In 



ttR^L 



R 7ri?4 ttR^L 



■"'t^i- 



(6) 



This equation shows that for fixed separation distance L in the limit cFstriT) -^ +0 the leading term is given by 
the surface tension contribution, while the next to leading term corresponds to the contribution of the Fisher 
topological term, whereas the second term on the right hand side of ^ is the smallest one. Therefore, it is 
evident that for small values of string tension (and large R) the main contribution to the total pressure and to 
its temperature derivative is given by the first term on the right hand side of (|6|. 

To calculate the total entropy density Stot of the cylinder let us parameterize the string tension as 

'y^trt" (7) 



where t = -^ > +0 and u = const > 0. From (7) it follows R 



go\n(L/Lo) 



and then for t — ^ the 
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entropy density Stot can be found from ^ and ^ as 

d Ptot\ ^ CTsurf 1 d dsurf 



Stot 



dT 



2 R Tnn t 



dominant since t^- 



ROT 



2Teo [go\n{L/Lo)\ t^-^l^ 



^ surf 



> 0. (8) 



This equation shows that at T = Tco the entropy density diverges for i^ < 2 and also that at the cross-over region 
the surface tension coefficient must be negative otherwise the system would be thermodynamically unstable since 
its entropy density would be negative. 
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Figure 2. Surface tension of normal paraffins as a function of temperature from the triple point to the 
critical point. The filled circles indicate the experimental data [30^. The lines are the different theoretical 
parameterizations [28^ . 

Clearly, the results (p!])-([8l) are valid for nonzero baryonic chemical potential /i up to the (tri) critical endpoint. 
The main modification in ([lf^([8| is an appearance of /i-dependences of p^(T, /i) and Tco{l~i) [6 . In the (tri) critical 
endpoint vicinity the behavior of p tot and Sfot is defined by the T-dependence of the surface tension coefficient. 

We stress that there is nothing wrong or unphysical with the negative values of surface tension coefficient, 
since (Jsurf '^ttRL in (IT]) is the surface free energy and, hence, as any free energy, it contains the energy 
part Esurf and the entropy part Ssurf multiplied by temperature T, i.e. Fgurf = Egurf — TSgurf tlHl flQ] . 
Therefore, at low temperatures the energy part dominates and the surface free energy is positive, whereas at 
high temperatures the number of configurations of a cylinder with large surface drastically increases and the 
surface free energy becomes negative since Ssurf > -^^-^- Moreover, the exactly solvable models with phase 
transition and cross-over [HI 13 [10] have region of negative surface tension coefficient and they clearly show that 
the only reason why the 1-st order deconfinement phase transition degenerates into a cross-over at low baryonic 
densities is the negative values of cFgurf at this region and the above results independently prove this fact. 

We believe that the same is true for many ordinary liquids otherwise one has to search for an alternative 
explanation for the disappearance of the 1-st order liquid-gas phase transition at the supercritical temperatures. 
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Of course, the experimental data in this region do not exists, but, nevertheless, there is indirect evidence 
for an existence of negative values of the surface tension coefficient at the supercritical temperatures. To 
demonstrate the validity of this statement we have to remind that the modern experimental data on the 
temperature dependence of the surface tension do not allow one to definitely conclude what is T-dependence at 
the vicinity of critical temperature Tc. In fact there are two alternative prescriptions [261127] 



Eotvos rule : 



Guggenheim rule : 



^ surf 

2 
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^ surf 
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p! 
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(9) 
(10) 



where pi is the temperature dependent particle density of the liquid phase. 
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Figure 3. The surface tension T (in N/m) in terms of cluster surfaces (first and third rows) and the surface 
tension V j p^' (in {N/m)/{mole/lY^'^) in terms of cluster number (second and fourth rows)) as a function of 
e = {Tc — T)/Tc for: quantum fluids (hydrogen and helium), noble gases (krypton and xenon) and more complex 
fluids (methane and water). The thin solid lines show data points ^ and the heavy dashed-dotted lines over 
the thin line show fits to F = Fq p/ e^^ and T / p^ = Fq £ according to Eqs. (10) and (N), respectively. 



After the Guggenheim work [26^ the prescription ( 10) became a dominant one [28 . Sometimes there appeared 
even confusions. Thus, in ^29| the authors determined the surface tension of water from the triple point to critical 
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point and parametrized it by the polynomial of 9-tli power (7surf{T) = ^ ai {Tc — T)\ but then the same 

authors refitted it to the prescription (10) [30]. Here in Fig. |2] which is taken from [28 we show the temperature 
dependence of some paraffins. As one can see for n-Pentane and n-Heptane the data on temperature dependence 



of surface tension near the critical point, indeed, may show the nonlinear behavior similar to (10), but for the 
n-Hexane and n-Octane one can see the linear T-dependence of Eq. ([9|! 

Therefore, in order to clarify this issue a few years ago a thorough analysis [3? of the high quality NIST 
data [33 was performed. Some of the results are shown in Fig. [s] which is taken from Ref. [32 . As one can see 
from Fig. Is] for most of the analyzed liquids the linear prescription ^ provides an essentially better fit with 
the only exceptions of xenon and methane. Therefore, our first conclusion is that for many liquids the rule 
([9| better describes the data than the rule (10). The second conclusion one can draw from this discussion is 
that naive extrapolation of the linear T-dependence ([9| of the surface tension coefficient cFgurf to supercritical 
temperatures T > Tc would lead to the negative values of the surface tension coefficient. Of course, one may 
think that asurf = for T > T^ like in the FDM [17 and SMM [22, 23, 24 , but in this case one has to explain 
the reason why the T-derivative of (Jsurf has a discontinuity at T = Tc while the pressure and all its first and 
second derivatives are continuous functions of it arguments in this region. 



4 The mysterious maximum of the lattice entropy and the gas of free tubes 

The considered configuration of the unbroken confining tube is only one of many other configurations accounted 
by the lattice QCD thermodynamics. However, in order to explain a mysterious maximum of the lattice 
entropy (see Fig [l]) it is sufficient to assume that the probability of the unbroken confining tube among other 
configurations measured by lattice QCD is W{L) ~ [L go\n{L/Lo)]~^^ i.e. in the limit T ^ oo it is negligible 
for any i^ 7^ 0. Then the contribution of the unbroken confining tube into the lattice free energy is small, since 
W{L)Fstr{L) ~ R~^ for t -^ +0 and R — > Riat — {Riat denotes the lattice size), but its contribution to the 
tube entropy 



W{L)Ss 



-W 



dFstr 

dT 



WL- 



-t^- 



WL- 



Tco L[^oln(T/To)]i- 



R 



2(1-^) 



R 



2(1-^) 



(11) 



is an increasing function of the tube radius Rioi v < 1. Clearly, if the available size of the lattice Riat would be 
infinite then the contribution of the unbroken tube would diverge, but for finite lattice size one should observe 
a fast increase at T ^ Tm- 



The physical origin of a singular behavior of the tube entropy ( 11 ) encoded in i^ < 1 is rooted in the formation 



of fractal surfaces of the confining tube in the cross-over temperature vicinity [6 . This can be clearly seen from 
the power ^ ~^^ of R on the right hand side of (11) which is fractal for any u ^ 2+^ where n = 1,2,3,... 
Moreover, the appearance of fractal structures at 1 = Tco can be easily understood within our model, if one 
recalls that only at this temperature the fractal surfaces can emerge at almost no energy costs due to almost 
zero total pressure ^. An explanation of the tube entropy decrease for t < is similar [6l[7]- It means that the 
fractal surfaces gradually disappear since for T > Tco the tube gradually occupies the whole available lattice 
volume. 

Here we also would like to consider a toy model based on the total pressure (IgI) of the confining tube, 
but for the non-static (or free) quark- ant iquark pair. In this case the parameter L should be considered as a 
free parameter which has to be determined from the maximum of the total pressure (p|. Finding from this 
condition the radius dependent separation distance Ly^{R) which corresponds to the most probable and the 
stable state of the free confining tube one has to substitute it into expression for the pressure ([6| and find the 
corresponding radius of the tube from the equation p = ptot{Lw{R)i R^ T) for the set of given external pressure 
p and temperature T. Clearly the determination of an extremum of the total pressure (|6| with respect to R 
with subsequent finding of the separation distance L should give the same result, but the first way is technically 
easier. 

Instead of the quantitative analysis of the resulting pressure ptot{Lw{R)i R^ T) which requires the knowledge 
of values of all constants, i.e. Tq, go^ t e.t.c, we prefer to discuss some qualitative properties of the model and 
show that it has two phases. Also we have to stress that such a model cannot be applied to high pressures 
(or high densities) directly because in this case the pressure of the system should account for the short range 
repulsion between the tubes. Therefore in what follows it is assumed that the gas of tubes has some low particle 
density p and, hence, one can consider this gas as an ideal gas with the pressure p = Tp. 

To determine the density p one has to maximize the pressure ptot first. From the vanishing derivative 
condition 



S ptot{L,R,T) 
S L 



90 



Tr 



ttR^L ttR^L^ 



In 



ttR'^L 



Tr 



ttR^L^ 



(12) 
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one can find the following equation for Ly^{R) 



-L/m — 



TrR^ 

90 



2 -In 



ttR^L^ 



(13) 



and show that it corresponds to a maximum of pressure. From Eq. (12) it is clearly seen that, in contrast to 



previous findings [HIT], the role of the Fisher topological term is a decisive one for an estabhshing Eq. (13). 
Consider a few limiting cases of Eq. (13). In the limit T ^ and finite R one gets 



L?. 



TrR^ 

90 



In 



ttR^tT 

Vo go e^ 



0, 



(14) 



which can be interpreted as a confinement of color charges. The same solution is true for the case of R ^ 



and finite T. In the other extreme T ^ oo (or R -^ oo) one obtains a different solution of Eq. (13) 



L°. 






(15) 



which again shows that for large R values the separation distance is small. This situation resembles what is 
observed in lattice QCD for the non-static color charges: the long tubes that connect such charges simply break 
up at some separation distances [Ij. 

From Eq. ( 13 ) one can show that the solution L^, is a monotonically increasing function of T, while for 
T ^ it always has a maximum as the function of R. Searching for the maximum of L^{R) ([13]) one can find 
the corresponding value of the radius Rmax and Lyj{Rjnax) 



-t^max 



goVpe 

tttT 

VoerT 
■^90 



(16) 
(17) 



which, evidently, obey the condition 7rR'^^^Luj{Rmax) = Vo e. The presence of the maximum of function Ly^{R) 
leads to an existence of two different radii for the same value of separation distance L, or in other words, there are 
two solutions of the equation L = L^{R). Clearly, the parameters of the maximum Rmax and L^{Rmax) given, 
respectively, by Eqs. (16) and (p!7[), separate the regions of these solutions. Evidently, the latter correspond 



to two phases of the gas of tubes which have different tube radius for the same separation distance L. The 
analysis of these solutions shows that there are many possibilities which strongly depend on the values of the 
involved constants Lq, go^ t and Vo^ whereas for low but non- vanishing temperatures one can show that the 
higher pressure corresponds to the phase of the tubes with smaller radius. This is clear because in case of low 
temperatures the leading contribution to the total pressure (pi) is given by the surface tension term ^^^^^^ 



R 



5 Conclusions 



In this work we discuss the most general relation between the tension of the color tube connecting the static 
quark- ant iqu ark pair and the surface tension of the corresponding cylindrical bag. Such a relation allows us to 
determine the surface tension of the QGP bags at zero temperature and, under the plausible assumptions that 
are typical for ordinary liquids, to estimate the temperature of vanishing surface tension coefficient of QGP bags 
at zero baryonic charge density as T^r = 152.9 ±4.5 MeV. Using the Fisher conjecture [17 and the exact results 
found for the temperature dependence of surface tension coefficient from the partition of surface deformations 
[18l[T9j[7j, we conclude that the same temperature range corresponds to the value of QCD (tri)critical endpoint 
temperature, i.e. Tcep = Tcr = 152.9 ± 4.5 MeV. Then requiring the positive values for the confining tube 
entropy density we demonstrate that at the cross-over region the surface tension coefficient of the QGP bags is 
unavoidably negative. Furthermore, analyzing the data on the temperature dependence of the surface tension 
coefficient of some ordinary liquids in the vicinity of the critical endpoint we conclude that the negative values 
of the surface tension coefficient of QGP bags are not unique, but they also should exist at the supercritical 
temperatures of usual liquids. We believe such a conclusion is worth to verify experimentally for ordinary 
liquids. 

Also we demonstrate that the long unbroken tube taken with a vanishing probability which generates a 
finite contribution into the lattice free energy may, under certain assumptions, provide a very fast increase of 
the lattice entropy of such configurations and, thus, it may explain the maximum of the tube entropy observed 
by lattice QCD. Additionally, we considered the non-static (free) tube and used the developed formalism to 
work out the model of the gas of free tubes. The performed analysis of such a model showed that there are two 
phases in this model which correspond to different radii of the tube for the same separation length L. Since this 
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toy model resembles some important features of the confining tubes observed in the lattice QCD, we believe 
it is important for QCD phenomenology and can be used to build up more elaborate statistical model which 
accounts for a realistic interaction between tubes. 
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